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Abstract. We present results on spectroscopic observa- 
tions of the fast-rotating active giant FK Comae, obtained 
mainly with the ESA-MUSICOS spectrograph at the Isaac 
Newton Telescope (INT) in 1996 and 1997 and also with 
the Aurelie spectrograph at the Observatoire de Haute 
Provence (OHP) in 1997. The profiles analysed are those 
of the Balmer, Ha and H/3, and Hei D3 (A 5876 A) fines. 

We analyse tfie Balmer line variability and phase be- 
haviour. We confirm these lines as highly variable, with 
excess emission that originates from extended structures 
and exhibits clear signs of rotational modulation. We have 
described the line profiles for two distinct states of activity, 
using different modelling approaches. Similar techniques 
were appfied to the He i D3 spectra. 

A large flare event lasting several days was detected 
in both Balmer lines and in the He i D3 line. The energy 
released during this flare in Ha is of the order of 10"^^ erg, 
making it the largest Ha fiare reported on a cool star. 

Our results confirm the extreme complexity of the cir- 
cumstellar environment of FK Comae. The activity level of 
this star is quite variable demanding different approaches 
to the line profile analysis. 

Key virords: Stars: activity- Stars: individual: FK Comae 
- Stars: late-type - Stars: flare - Line: profiles - circum- 
stellar matter 



1. Introduction 

FK Comae Berenices (HD 117555), with its extreme ac- 
tivity, is a key target to study in the solar-stellar con- 
nection context. It is a rapidly rotating and apparently 
single G5 II giant, even though the spectral classification 
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can go as late as G8 HI, due to the fast rotation rate 
and spectral oddities. This star has a rotational veloc- 
ity of V sin i = 162.5 ± 3.5 km s^^ and radial velocity 
-24 ± 3 km s"^ (Huenemoerder et al. 1993 hereafter ref- 
ered to as HRBN) . The photometric and rotational period 



is 2.400 day (Jetsu et al. 1993D . This gives a projected stel- 
lar radius of 7.7 ± 0.2 Rq (HRBN). As a member of the 
HR 1614 old disk aggregate, its distance and mass are ap- 
proximately known: d ^ 200-300 pc and M ^1.5 Mq 



(Eggen& Ibcn |l989D . 

Some spectral peculiarities of FK Com were first ob- 



served by Merrill (1948). Amongst the oddities present in 
the spectra of this star can be included a strong, broad 
and erratically variable Ha emission (e.g. Ramsey et al. 
1981 ), strong UV chromospheric emission lines reminis- 



cen t of t hose observed in RS CVn binaries (Bopp & Sten- 
cel 1981 ), a quasi-sinusoidal light curve (Holtzman &Na- 
tions 1984 ) and high X-ray luminosity (Walter 1981). 
Some of these spectral characteristics led Bopp & Stencel 
( 1981 ) to define a new class of variable stars with FK Co- 
mae as the prototype. The FK Comae-type variables were 
thus defined as rapidly rotating, single, G-K type giants 
with strong chromospheric and transition region emission. 
Stars belonging to this class are extremely rare. The space 
density of these stars was estimated by Collier ( |1982D as 
2 xl0~* pc~'^, indicating, when compared with the space 
density of the rare Hertzsprung-gap G giants, that the 
evolution through this part of the H-R diagram must be 
extremely rapid. 

Typical values of rotation rate from G-type stars are 
< 10 km s^^ for single stars and 30-50 km s~^ for stars 
in close binary systems. The extreme rotational velocity, 
near breakup, of FK Comae raises interesting questions on 
the possible evolutionary status of this star. Two scenarios 
were advanced to explain the observed rotational and ac- 
tivity characteristics: accretion from a very small unseen 
companion or product of the recent coalescence of a close 
binary system. Based on radial velocity variations mea- 
surements, McCarthy & Ramsey ( 1984 ) proved the first 
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hypothesis very unhkely, by imposing severe Umits to the 
mass of this unseen companion. The currently accepted 
theoretical scenario is that FK Comae-type stars are re- 
cently coalesced binary systems. According to Webbink 
( |1976D , W UMa systems, due to the thermal instabilities 
caused by the primary's departure from Main-Sequence, 
evolve into configurations that can only be stabilized by 
a big loss of angular momentum, coalescing into a single 
star in a way that is mostly independent of their evolution. 
The rapid rotation would arise from the angular momen- 
tum of the former binary and the activity level from a 
rotation-driven magnetic dynamo (Ramsey et al. 1981 ). 



The descendents of these binaries would be ind istinguish- 
able from single stars. Welty & Ramsey ( 1994 ) presented 
evidence that conforms with such scenario, based on the 
analysis of the width and shape of the spectral lines. 

In FK Comae, the most striking spectral feature 
is the Ha emission line. This line profile is extremely 
broad, asymmetric and presents variations on different 
time scales. The same characteristics are present in higher 
Balmer lines, even though appearing less extreme due to 
the different optical thickness of the lines. The phase mod- 
ulation in these lines indicates that the emitting extended 
material is corotating with the star. Based on the Ha to 
H/3 ffux ratios, HRBN have proposed that the emission 
arises in structures similar to solar prominences. 

We present here a model of the Ha and H/3 line profiles 
for a quiet phase of the circumstellar environment, derived 
from a data set obtained in May 96. The model includes 
a 3D gaussian spatial distribution of emitting matter, ex- 
tending to several stellar radii, and an optically thick near- 
stellar shell of absorbing material. The data sets from May 
and June 97 represent a more active state, demanding a 
different modelling technique. The H/5 line profiles present 
a complex structure that we model as the superposition of 
gaussian-shaped components. The same components were 
found to be present in the Ha line profiles, superposed 
on a more distributed "bulk" emission. We were able to 
track these components in the velocity vs phase diagram 
for approximately 3 rotational periods. A large fiare event 
was detected in our May 97 data set, in the Balmer and 
He I D3 lines, and we computed the corresponding energy 
budget. 

2. Observations and data reduction 

The MUhi-SIte Continuous Spectroscopy (MUSICOS) 
project was developed with the aim of providing contin- 
uous spectroscopic coverage (in a large wavelength range 
and with instruments as similar as possible) for a wide 
range of stellar objects (e.g. Foing & Catala [l991 : Foing et 
al. |1994bt Catala et al. |1990|, |1993|; Huang etal.|1995D. Fol- 



lowing these orientations, the first fiber-fed echelle cross- 
dispersed MUSICOS spectrograph was designed and built 
in Meudon-Paris Observatory (Baudrand & Bohm 1992 ) 
and a replica was built at the ESA Space Science De- 



partment in ESTEC and adapted for the INT. The ESA- 
MUSICOS spectrograph was installed at the 2.5 m Isaac 
Newton Telescope (INT), at the Observatorio del Roque 
de los Muchachos, La Palma, Spain. 

The commissioning campaign of this instrument took 
place from April 29th to May 6th. Our first data set was 
obtained in this campaign. The resolving power achieved 
was A/AA ^ 35 000. The second data set was obtained 
in May 97 also at the INT with the same instrument. 
The CCD detector used in both runs was a Tektronix, 
1024 X 1024 pixels, with a readout noise of ~ 5 e~. To 
reduce the INT data, we used the MIDAS Nov. 94 and 
Nov. 96 echelle reduction packages. The reduction pro- 
cedure used was standard for echelle spectra, with some 
problems in the correction for the blaze function due to 
the low level of counts of the fiat-field exposures in the 
bluer part of the spectra. 

In addition to these INT observations, we obtained in 
June 97 more Ha spectra at the Observatoire de Haute 
Provence (OHP), in France. We used the Aurelie spectro- 
graph at the 1.52 m coude telescope, with resolving power 
of 30 000. The reduction of this data set was done using 
the on-line reduction facilities available on site. 

In Table 1 we give the summary information on 
all the spectra obtained. The ephemeris used for 
the phase computations is from Chugainov (1976): 
2 442192.345 + 2.400E. 

2.1. Artificial template for excess spectra 

As we are mainly interested in the study of phenomena 
related with the excess activity, one more standard re- 
duction step was applied to all the spectra. As long as 
the line profile has the same shape over the stellar disc, 
a rotationally broadened ffux profile can be obtained by 
convolving the ffux profile of a non-rotating template star 
with the rotational profile of the observed star. As a tem- 
plate, we used k Hercules (HD 145001), a G8 III type star, 
still adequate for differential comparison. We have artifi- 
cially broadened the spectrum of this low activity star 
to simulate the rotational broadening of FK Com. This 
reproduces adequately the photospheric lines and serves 
as baseline for the chromospheric lines. After normaliza- 
tion, these model spectra were subtracted from the ones 
of FK Com, for the spectral orders of interest. 

3. Results from May 96 line profiles 

3.1. Low activity Balmer extended emission 

From 1st May 96, we describe 4 FK Com spectra of Ha 
and H/3. These spectra are symmetric and the measured 
equivalent widths (respectively -2.4 and 1.0 A) are low 
when compared to other available measurements (both 
from the literature and from our measurements). This sug- 
gests that these spectra represent a low activity state for 
this star. 
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Table 1. List of 1996 and 1997 FK Comae observations 



HJD* 


Date 




UT* 


tr 


site 


0'" 


2 450 205.38 


96 May 


1 


21:08:20 


1800 


INT 


.764 


2 450 205.41 


96 May 


1 


21:45:40 


1800 


INT 


.775 


2 450 205.44 


96 May 


1 


22:38:41 


1800 


INT 


.790 


2 450 205.46 


96 May 


1 


23:09:12 


1800 


INT 


.800 


2 450 583.47 


97 May 


14 


23:20:57 


900 


INT 


0.343 


2 450 583.48 


97 May 


14 


23:39:03 


1800 


INT 


0.345 


2 450 584.36 


97 May 


15 


20:37:44 


1500 


INT 


0.712 


2 450 584.50 


97 May 


16 


00.05:44 


800 


INT 


0.772 


2 450 585.37 


97 May 


16 


20:49:30 


900 


INT 


1.132 


2 450 585.39 


97 May 


16 


21:06:29 


1200 


INT 


1.137 


2 450 585.39 


97 May 


16 


21:27:55 


1200 


INT 


1.143 


2 450 585.41 


97 May 


16 


21:48:20 


1200 


INT 


1.149 


2 450 587.38 


97 May 


18 


21:08:37 


1200 


INT 


1.971 


2 450 587.40 


97 May 


18 


21:29:41 


1200 


INT 


1.977 


2 450 588.38 


97 May 


19 


21:06:12 


1200 


INT 


2.387 


2 450 588.39 


97 May 


19 


21:26:56 


1200 


INT 


2.393 


2 450 589.38 


97 May 


20 


21:08:28 


1200 


INT 


2.804 


2 450 589.39 


97 May 


20 


21:29:14 


1800 


INT 


2.810 


2 450 590.38 


97 May 


21 


21:08:30 


1800 


INT 


3.221 


2 450 590.40 


97 May 


21 


21:38:54 


1200 


INT 


3.230 


2 450 617.38 


97 June 17 


21:04:56 


1200 


OHP 


14.470 


2 450 617.43 


97 June 17 


22:12:53 


1200 


OHP 


14.489 


2 450 617.47 


97 June 17 


23:16:24 


1200 


OHP 


14.508 


2 450 623.39 


97 June 23 


21:20:55 


1200 


OHP 


16.974 


2 450 623.41 


97 June 23 


21:44:47 


1200 


OHP 


16.981 


2 450 623.42 


97 June 23 


22:08:28 


1200 


OHP 


16.988 


2 450 623.44 


97 June 23 


22:33:26 


1200 


OHP 


16.996 


2 450 623.46 


97 June 23 


22:58:06 


1200 


OHP 


17.003 


2 450 623.48 


97 June 23 


23:22:20 


1200 


OHP 


17.010 



* Times given at the start of the exposures. 

** Exposure time in seconds. 

*** A rotation-number ordinate was added to the 4> value of 

the spectra of 1997 to represent the diflerent rotation periods 

covered. 



The difference spectra show emission extending up to 
~ 780 km s~^ in Ha and ~ 450 km s~^ in H/3. We found 
that in both lines the wings of the extended emission, out- 
side V sin i, can be fitted by gaussian profiles, as it can be 
seen in Fig. n^. The full width at half maximum of the fitted 
gaussian profiles are 642 ± 12 km s~^ and 381 ± 20 km s^^ 
respectively for Ha and H/3. Assuming that the velocities 
correspond to corotating material, we derive the maxi- 
mum extent of the emission, ^ 4.75 iZ, sin i for Ha and 
'^ 2.66 i?, sin i for H/3 (i?* sin i is the projected stellar 
radius). 

Using the Cousins-Bessel magnitudes for FK Com, 
TaR = 7.66 and m^ = 9.07 (Hohzman & Nations 1984) 
and the corresponding flux cahbrations (Bessel 1979) the 
Balmer line fluxes were computed. The ratio between the 
fluxes of the Ha and H/3 for the fltted gaussian is at max- 
imum emission ~ 4, and '^ 7 for the integrated flux, rem- 
iniscent of solar-like active prominences (Oliveira et al. 
1998a). However, as we derived different widths for Ha 





and H/3 indicating a different extension of the emitting 



Fig. 1. The Ha (top) and H/3 (bottom) average excess 
spectra on 1st May 96, as well as the gaussian profiles 
(dotted line) , fitting the wings of the extended emission. 
The obscuration by the stellar disc and the absorption of 
the near-stellar shell are displayed (dashed line) and also 
the effect of the extra absorption of the background stellar 
profile (full line). 



plasmas, these global ratios may have little direct physi- 
cal meaning. 

3.2. Modelling of circumstellar emission 

We have modelled the emission wings by assuming an 
effectively thin, corotating emission and a 3D gaussian 
distribution of the source function. Under such assump- 
tions, the measured full width relates directly to the scale 
height of the emission. We have performed numerical in- 
tegration of such a 3D distribution at each velocity. This 
characterization represents well the observed wings, but 
differs strongly in the central absorption profile, within 
^ ± 1.3 u sin i from the rotational axis. 

The central absorption, superposed to the emission, 
can partly be interpreted as the presence of the star itself 
and the effect of obscuration from the stellar disc. The 
analysis of the profiles shows that this absorption is in 
fact slightly broader than 2 x v sin «, in the two Balmer 
lines, suggesting the presence of a near-stellar shell of ab- 
sorbing material. The radius of this absorption shell Rabs 
is one of the parameters to be adjusted in this modelling 
approach. The presence of this shell also contributes with 
another extra absorption, as it absorbs the background 
stellar profile by a factor e. 

The numerical integration of the 3D gaussian emission 
is computed at a given velocity, removing the volume ob- 
scured by the star and the absorption shell. The modelled 
profile obtained by this method is then corrected for the 
extra absorption, by subtracting the function e x I^, where 
If, is the background stellar spectrum. 



Oliveira & Foing: Circumstellar emission and flares on FK Comae 



The resulting profile is compared with the difference 
spectrum in Fig.^ The values for the parameters that 
give the best agreement with the observed line profiles are 
R fj^« = 1.25 i?, sin i, Rf^f = 1.1 R^ sin « and £ = 0.225. 

The fractional depression of the line profiles is approxi- 
mately the same in these two Balmer lines, indicating pos- 
sibly an optically thick absorption. Under this assumption, 
the value of e can give us an estimation of the fraction of 
the stellar disc obscured by the near-stellar shell; accord- 
ing to this model we have a disc coverage of about 23 %. 

The residuals, after this axisymmetric model contribu- 
tion is removed, have an Ha to H/3 flux ratio of the order of 
4, indicating an optically thicker source than the extended 
emission. Their projected velocities are within ± v sin i. 
Thus, they can be interpreted as plage-like structures and 
they account for ^^ 10 % of the total emission. 

4. Results from May/June 97 line profiles 

4-.1. Multi- component analysis 

In May 97, FK Comae spectra of Ha, H/3 and Hei D3 
were obtained at the INT. This data set is shown in the 
4 upper rows of graphs in Fig. g. Also in Fig. g, in the 
lower row we show the Ha spectra obtained at the OHP, 
3 in June 17th and 6 in June 23rd. Each pannel shows the 
measured spectral lines for representative times, after re- 
moving the rotationally broadened template spectra (each 
spectrum is an average of the spectra taken in each night). 
We have tried to model these spectral lines in a similar 
way as it was done for the May 96 data set. Clearly, that 
axisymmetric model does not apply to these May/ June 97 
Balmer line profiles, very asymmetric, variable and with 
a more complex component structure. Thus, we used a 
multi-component approach, that consists in fitting a vari- 
able number of discrete gaussian components to each of 
the spectral lines, representing emitting structures. We 
also analysed the phase behaviour of these structures. We 
plot these components in velocity vs phase ((/)) . We fit the 
velocity curves that describe the phase behaviour of the 
emitting structures as they corotate with the star, by an 
equation of the type: 



R 
R» 



X v^q sm I cos a sm( 



Vc X sm 



where Veq is the equatorial rotational velocity, 9 is the 
stellar latitude, (pQ the stellar longitude and Vc is the ve- 
locity excursion in km s^^. The two Balmer lines were 
analysed in parallel to investigate their similarities. The 
same multi-component approach was also used for the He i 
D3 line. 



^.2. The Balmer lines fitted components 

Figure |^ shows an example of how the profiles were fit- 
ted with gaussian emission components. Comparing the 



two Balmer line profiles in Figs. and it is clear that 
the same components seem to be present in both spectra. 
This becomes even more evident when the average Ha 
and H/? profiles are subtracted from the individual pro- 
files. The analysis of these average profiles hints that the 
H/3 emission is somehow dominated by discrete compo- 
nents, while in Ha these same components are masked by 
a more stationary, distributed, "bulk" emission. 

Ha and Hf? fitted components 
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Fig. 3. An example of the multi-component fit to the 
Balmer line profiles (from top to bottom H/? and Ha). 
It is clear that the same components are present in both 
spectral lines, except for the higher velocity ones, that are 
masked in the noise in H/3. 

In Table 2 we give the parameters for each of the gaus- 
sian components fitted to the Balmer line profiles. Ha and 
H/3. In Fig. 0, we show several plots that describe their 
behaviour: from top to bottom, the peak intensity vs ve- 
locity, and velocity vs phase for both lines. In order to 
highlight their phase behaviour, we also plotted (in the 
two lower graphs) relevant velocity curves. Figure shows 
the velocity evolution of the components in time. We have 
estimated a maximum error in the velocity determinations 
of about 50 km s^^. The plots for these two Balmer lines 
were analysed together in order to clarify the structures' 
position and evolution. 
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Date 97 





Vq/ii sin i 


(FWHM)„/^ sin I 


A, 


Vfj/v sin I 


(FWHM)^/^ sin I 


A^ 


identification 


14 May 


0.343 


-3.3 


1.0 


0.08 


+ 


* 


+ 






0.343 


-1.2 


2.4 


0.35 


-1.1 


1.8 


0.13 


A 




0.343 


1.5 


3.2 


0.34 


1.8 


1.7 


0.06 


B 


15 May 


0.712 


-1.8 


2.6 


0.40 


-1.6 


1.6 


0.06 


B 




0.712 


-0.1 


1.0 


0.19 


0.0 


0.6 


0.10 






0.712 


1.1 


2.1 


0.40 


1.2 


1.2 


0.08 


A 




0.712 


3.4 


1.1 


0.07 


+ 


* 


* 




15 May 


0.772 


-2.2 


2.2 


0.29 


-1.7 


1.4 


0.06 


B 




0.772 


0.6 


3.3 


0.53 


0.9 


1.8 


0.16 


A 


16 May 


0.132 


-1.9 


2.9 


0.67 


-1.5 


2.5 


0.21 


A 




0.132 


-0.3 


1.2 


0.31 


-0.4 


0.6 


0.11 






0.132 


1.4 


3.0 


0.47 


1.9 


2.5 


0.09 


B 


18 May 


0.971 


-0.6 


2.9 


0.32 


-0.9 


1.6 


0.06 


C 




0.971 


1.5 


1.3 


0.25 


1.4 


1.0 


0.10 






0.971 


2.3 


2.9 


0.10 


2.9 


0.3 


0.04 




19 May 


0.387 


-3.4 


0.7 


0.07 


* 


* 


* 






0.387 


-0.9 


2.2 


0.29 


-0.8 


2.1 


0.12 






0.387 


0.3 


1.0 


0.37 


0.2 


0.6 


0.14 


A 




0.387 


1.5 


1.2 


0.37 


1.3 


1.5 


0.14 


B 




0.387 


2.9 


1.3 


0.09 


* 


* 


* 




20 May 


0.804 


-2.8 


1.9 


0.11 


-2.8 


1.6 


0.02 


C 




0.804 


-1.5 


0.9 


0.22 


* 


* 


* 


B 




0.804 


0.4 


1.8 


0.50 


1.0 


1.3 


0.08 


A 




0.804 


1.9 


3.4 


0.12 


+ 


* 


* 




21 May 


0.221 


-1.7 


2.6 


0.26 


-0.7 


1.2 


0.08 






0.221 


0.9 


2.5 


0.38 


1.6 


1.0 


0.07 


B 




0.221 


3.1 


1.3 


0.07 


2.9 


1.0 


0.03 


C 


17 June 


0.470 


-3.4 


0.4 


0.02 


** 


** 


** 






0.470 


-1.1 


3.5 


0.43 


** 


** 


** 






0.470 


-0.2 


0.7 


0.16 


** 


** 


** 






0.470 


1.3 


1.4 


0.20 


** 


** 


** 






0.470 


2.9 


1.3 


0.07 


** 


** 


** 




23 June 


0.988 


-1.7 


2.8 


0.27 


** 


** 


** 






0.988 


0.2 


0.5 


0.04 


** 


** 


** 






0.988 


1.6 


3.0 


0.39 


** 


*+ 


** 





Table 2. Parameters for the fitted Ha and H/3 components for the INT May 97 and OHP June 97 campaigns. Vq 
and (FWHM)q are respectively the velocity position of the maximum and the full width at half maximum (in units 
of V sin i = 162.5 km s~^) for each of the fitted gaussian components in the Ha excess spectra. The velocity errors 
are typically 0.1-0.3 v sin i. Aq is the amplitude (in continuum units) for each of the components, v^g, (FWHM)^ 
and A^ are the corresponding quantities for the H/3 components. Some Ha components do not have measurable H/3 
counterparts (indicated by *). For the June 97 data set (**), no H/3 observations were made. The capitals A, B and 
C in the last column identify some of the components which behaviour is described by the velocity curves plotted in 
Figs. 4 and 5. 



4.2.1. Velocity curves analysis 

As it can be seen in Fig. ^, three structures could be fol- 
lowed in phase, corresponding to va = 245 x sm{(j)+0.62), 
vb — 285 X sm{ip) and vc — 487 x sin((/)), both in Ha and 
H/3. It was observed that most of the components were de- 
tected within approximately ± 2 w sin i from the rotational 
axis, in agreement with Welty et al. (1992) that placed the 
emitting material at < 2 i?,. 

Similar component analysis was carried out with the 
Ha spectra obtained at the OHP in June 97. Two repre- 
sentative spectra were fitted and the velocity positions of 



the components were compared with what was obtained 
for May 97. Even though some profiles appear similar 
(Fig. |, for phase 1.971 and 16.974), the previously fit- 
ted velocity curves do not seem to describe the phase 
behaviour now observed (Fig. middle graph). In par- 
ticular, one can notice the appearance of a new compo- 
nent near-zero velocity at phase ~ 0.47, compared with 
one month earlier. However, the poor phase coverage in 
June 97 makes a detailed comparison ambiguous. We also 
find an increase of Ha equivalent width of 2 A in June 97 
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compared to the quiet level of May 97, for the two phases 
shown. 



4.2.2. Comparison of Ha and H/3 components 

Most of the components detected in Ha have their coun- 
terparts in H/3. Still, some higher velocity Ha components 
are not seen in H/3. This is due to the better detection sen- 
sitivity in Ha as the signal-to-noise ratio in H/3 is lower, 
but also to the presence of the "bulk" emission in Ha that 
extends to higher velocities (as it was also seen in the quiet 
state in May 96 we modelled earlier). 

A slight mismatch can be observed in the velocity po- 
sition of the components in both lines. Intrinsic measure- 
ment uncertainties and physical phenomena combine to 
create this effect. The origin of the emission is obviously 
not in point sources, their vertical extent making difficult 
to determine the center of mass of the matter distribution, 
due to projection effects and emission gradients in the ex- 
tended structures. This lack of exact agreement between 
the two lines is then justified given the simplified geomet- 
rical model adopted and the ignoring of radiative transfer 
effects. 

4.3. The He I D3 line profiles 

The He I D3 (A 5876 A) line is very weak or absent in 
inactive stars' spectra. The presence of this line in cool 
stars like FK Com indicates the existence of non-radiative 
heating processes, since the line cannot be excited at pho- 
tospheric temperatures. In FK Comae, this line has two 
components: an excess absorption and near-surface emis- 
sion wings. A non-uniform shell of excited He I could ex- 
plain this line profile, with the corona as source of ionizing 
radiation as it is observed in the Sun and active stars. 

A similar multi-component study was made for the He i 
D3 line. All the spectra present an absorption at near-zero 
projected velocity, even though its strength is variable. 
Emission wings were detected above the limb, confirming 
the existence of a slightly extended shell of excited He i. 
The variability observed in the spectra can represent a 
back-radiation of a nonuniform coronal emission onto the 
upper chromosphere. The modelling technique we used did 
not allow us to take advantage of the rotational phase 
information that might be present in the He I data set, 
due to the complexity of the profiles. 

4.4- The 16 May 97 long-duration giant flare event 

A large flare event was detected in these two Balmer lines 
and also in Hci D3. The flaring material seems to extend 
from the stellar surface up to at least 3 stellar radii. 

The two spectra where we detected an increase in the 
excess emission, what we call a pre-flare spectrum and 
a flare spectrum, (respectively UT date 16.0 and 16.8 in 
Fig. II) are marked by "i^' in Figs. ^ and ^. 
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Fig. 6. The equivalent width (EW) variations for Ha and 
H/3, for May/June 97, both in UT date and phase. The 
lower values in each plot represent H/3. The dotted line 
indicates the lower Ha EW value measured for the aver- 
age spectrum from May 96. The symbols * represent the 
measurements from the OHP. The errorbars are 0.1 and 
0.2 respectively for the Ha and H/3 measurements. 



If these spectra represent the same flaring structure 
viewed at different phases, we can predict a velocity 
curve allowing to identify the post-flare structures at later 
phases. Such a velocity curve is also plotted in Fig. ^, cor- 
responding to Vc ^ 245 kms~^ (R = 1.50 i?* sin i) and 
00 = 0.612. At five other phase positions the same struc- 
ture is identified. Indeed, this structure was first detected 
on May 15.8 and it was still present on May 21.8, when 
our run ended. When transiting in front of the stellar disk, 
near zero projected velocity, this structure (extended ver- 
tically) would have smaller projected dimension, and in 
the spectrum we should detect a narrow peaked compo- 
nent. This is in fact observed in the spectrum correspond- 
ing to UT date 19.8 in Fig. |. 

This flare event detected in the Balmer lines caused 
an increase in equivalent width of 7 ± 1 A in Ha and 
1.2 ±0.2 A in H/3 (Fig. ||). In the case of H/3 the measure- 
ments were very difficult, the flare increase in equivalent 
width being the only variation detected above the error- 
bars. Using the Cousins-Bessel magnitudes for FK Com 
and assuming a distance of 215 pc (HRBN) some flare 
energetics were computed. The calculated energy released 
at flare maximum is 8 ± 1 xlO'^'^ erg s~^ for Ha and 
1.1 ± 0.1 xlO^i erg s~^ for H/3. The Ha to H/3 flux ra- 
tio is about 7, perhaps indicating much lower densities 
when compared with solar flares. Again the interpretation 
of such global ratios is dubious in this type of analysis. 

The total energy released during the flare for Ha is 
11 ± 3 xlO^^ erg. The uncertainties in the intrinsic flare 
radiative losses for Ha come mainly from the choice of 



Oliveira & Foing: Circumstellar emission and flares on FK Comae 



reference quiet spectrum and residual emission from ac- 
tive regions around the flare area. Comparing this en- 
ergy with the one reported on the HR 1099 giant Ha 
flare, ~ 4 xlO^^ erg (Foing et al. 1994a), this seems 



to be possibly the largest Ha flare reported on a cool 
star. What is also remarkable is the long duration of this 
flare event (total equivalent flare duration approximately 
2.5 ± 0.5 xlO^ s). 

The ratios of the continuum UBVRf counterparts to 
Ha, for similar flare events, are reported to be from 120 



(Avrett et al. |1984| ; Mac hado et al. |l980t Houdebine [1992D 
up to 340 (Foing et al. |1994aD. The continuum UBVRI 



counterpart of the FK Comae flare could therefore amount 
to losses of 1.3 xlO^^ erg (for a flare ratio UBVRI/Ha of 
120). This can be compared to total radiative budgets 
for the largest report ed cool star flares, 1.2 xlO^^ erg on 



HR 1099 (Foing e t al. |l994aD or 1.8 xlO^^ erg on YY Men 
(Cutispoto et al. 1992 ). Thus, we can speculate that we 
may have observed an exceptionally energetic flare, but we 
have no independent photometric measurements to quan- 
tify the continuum counterpart. 

An enhancement of He i emission was measured at pre- 
flare and flare maximum with equivalent width 0.08 and 
0.19 A respectively and at a velocity consistent with Ha 
and H/3 components. It was not measurable during the 
gradual phase, and in particular not in the post-flare spec- 
tra of May 19.8. This is similar with He I measurements 
in solar and red dwarf flares, which show a more impul- 
sive behaviour of this line, while the Balmer lines show 
enhancements both in impulsive and the gradual phase 
(Foing 1989). It seems also that the Hei D3 flare maxi- 
mum was reached before the Ha flare maximum. 

The changes observed in the large flare component 
width might also be associated with large kinetics (turbu- 
lence and flows) related with the flare event (Foing 1989| , 
Foing et al. |l994aD. 



5. Conclusions 

Once again, the analysis of the line profiles gives us abun- 
dant evidence on the complexity of the circumstellar en- 
vironment of the peculiar giant star FK Comae. 

We have modelled empirically the Balmer line profiles 
for this star. For a lower activity level, the model pictures 
a 3D gaussian emission to which a near-stellar thick ab- 
sorption shell is superposed. 

For the more complex structure of the 1997 profiles 
a different modelling technique is applied. We have suc- 
cessfully identified and described the phase behaviour of 
emission components. The H/3 spectral line was described 
as the superposition of several (gaussian) emission com- 
ponents. The velocity vs phase diagram allows to trace 
the same emitting structures as they corotate with the 
star. In the case of Ha, the profile is interpreted as a 
quasi-stationary distributed "bulk" emission, to which dis- 
crete components, lasting several days, are superposed. 



The same components were detected in both Balmer line 
profiles. 

A strong flare event was detected in these two Balmer 
lines and also in Hei D3. The flaring material seems to 
extend from the stellar surface up to at least 3 stellar 
radii. The evolution of the flaring structure was described 
by a velocity curve corresponding to Vc '^ 245 kms~^ or 
R = 1.50 R^, sin i from the rotational axis. The total 
energy released during the flare in Ha is ^ 10^^ erg, the 
largest Ha flare energy reported on a cool star. 

The data set of June 97, as it presents longer time 
intervals continuously covered, provides some more infor- 
mation. Short time-scale variations were observed, that 
cannot be associated with rotational modulation. Several 
variable features are visible; they can be interpreted as 
emission bursts with a time scale of about 1 hour that can 
be related to the excitation and recombination of giant 
magnetic loop structures around the star (Oliveira et al. 



1998b) 



The He i D3 line is always present in the spectra, with 
absorption and emission wings components, as a clear in- 
dicator of the activity in this star. An absorption near-zero 
projected velocity is always present with variable inten- 
sity. The profiles also present near-stellar emission wings, 
particularly enhanced during the described flare event. 

The need is clear for a multi-site continuous data 
set to disentangle the intrinsic variations from rota- 
tional modulation. We have organized such a campaign 
in February/March 1998, involving the following tele- 
scopes/observatories: INT (with the ESA-MUSICOS spec- 
trograph), Telescope Bernard Lyot, France (with the 
MUSICOS spectropolarimeter), David Dunlop Observa- 
tory (Canada), Russian Special Astrophysical Observa- 
tory (SAO) and McDonald Observatory. The phase cover- 
age thus obtained represents a major improvement in the 
analysis of the circumstellar emission on FK Comae. This 
will be reported in a subsequent paper (Oliveira et al. in 
preparation) . 
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Fig. 2. FK Comae excess spectra from May 97 at INT (upper four rows in the graph) and June 97 at OHP (bottom 
row). In the upper graphs (INT), we plot the Ha, H/3 and Hei D3 spectra. The Hei D3 spectra were multiphed by a 
factor 2, to improve features visibiUty. The Ha and H/3 spectra were displaced vertically by 0.3 and 0.15 respectively 
for clarity. In the lower row we show the Ha spectra obtained at OHP, using the sa me plo tting scale. In each graph 
the phase and UT date are also indicated. The ephemcris used is from Chugainov ( 1976 ): 2 442192.345 + 2.400E. 
The integer rotation-number ordinate describes the continuity of the data set, with the first spectrum of 14 May as 
reference. Velocity components are noticeable in the spectra and their variations in velocity and intensity were followed. 
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Fig. 4. Phase behaviour of Balmer hne components. From top to bottom: peak intensity plotted against phase with 
the size of the symbols representing the corresponding intensity, • represents Ha and o H/3 from INT May 97 and o Ha 
from the OHP June 97; velocity evolution of the Ha components with phase; velocity evolution of the H/3 components 
with phase. In the bottom two graphs, three velocity curves (A, B and C) were overplotted, corresponding to different 
distances from the rotational axis and longitudes: va = 245 x sin((/) + 0.62), v^ — 285 x sin((/)) and vc — 487 x sm{<j>). 
The dotted part of these curves represent the occultation of the structures by the star. The points indicated with "F' 
are flare spectra. The dashed lines represent ± v sin i from the rotational axis. 
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Fig. 5. Time behaviour of Balmer line components in May 97 (same plotting code as Fig. }^. The velocity curves of 
the previous figure were also overplotted. 



